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Pressure-Induced Phase Transition to a Novel Spin State in Striped Nickelates 
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We analyze pressure effects on stripe states within a selfconsistent Hartree-Fock calculation 
for a model of a striped nickelates. The results show a transition induced by high pressure and 
predict possible new spin states. We describe characteristics in the phonon excitations at the 
predicted transition, based on a real-space random phase approximation. 
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Much recent attention has been devoted to stripe 
patterns, 1-5 ' which are a typical spin and charge or- 
dering in many transition metal oxide compounds: 
La 2 -xSr x Cu0 4 , 6 ' 7) La2-xSr x NiO 4 , 8 " 10) etc. The stripe 
state has been intensely studied due to possible implica- 
tions for high-temperature superconductivity. 

To investigate the stripe state, many kinds of exper- 
iments have been carried out under a variety of condi- 
tions: different hole dopings, several kinds of impurities, 
etc. Neutron scattering is a powerful tool to detect stripes 
and such experiments have revealed important properties 
of the stripe state. For example, hole-doping effects on 
stripes has been studied by neutron scattering, in terms 
of the doping dependence of the stripe periodicity and 
orientation. 11 ) 

In addition to the doping effects, pressure should 
yield important information on the stripe states. How- 
ever, a pressure study is not readily compatible with 
neutron scattering experiments. Due to the difficulties of 
a high-pressure study on striped compounds, not many 
experiments have been performed. Nevertheless, high- 
pressure effects open up possibilities for finding new 
features and for understanding the properties of stripe 
states. Therefore, theoretical work is desirable to give 
guidance for high-pressure experiments. Here, we pro- 
vide a prediction of a new stripe state induced by high 
pressure based on a Hartree-Fock (HF) calculation for a 
model of striped nickelates with i-hole doping. 

We consider the hamiltonian 



H =TCo + Wish + Ho, 
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where i and j indicate a unit comprised of a Ni and two 
O sites, and m and n denote the orbitals (A = d x 2_ y 2, 
B = c? 3x 2„ r 2, p x or p y ) or specify an O ion (O x or 
O y ) when we consider ionic motion. Note that the sum- 
mations in 7^Ni (Ho) are taken only on the Ni (O) 
sites. The Ni on-site energy e, is equal to —A ± E z , 
where A = e p — and E z is the energy splitting of 
A and B orbitals caused by the crystal field. The term 
J' is a pair hopping, and its energy J' is set equal to 
the Hund exchange energy J for simplicity. tij^ mn (uk) 
is non-zero only for nearest neighbor hoppings and de- 
pends on the displacement of the involved oxygen-ion 
Uk through the electron-phonon coupling constant a: 
auk for tij, mn (u) ^ cases. For sim- 
plicity, we restrict the displacements of O ions only in 
the Ni-0 bond direction. 

We solve the model selfconsistently within a HF ap- 
proximation in a supercell of size 6x6 with periodic 
boundary conditions. To simulate the pressure change, 
we introduce a pressure parameter A, which is set as A = 
1.0 for ambient pressure and increased for higher pres- 
sure. The hopping energy is varied through A. We take 
the following parameter set: 12 ' 13 ) tA P — A, tB P = 
A = 9, Ez = 1, J = J' = 1, and U = 4 in energy unit 
i eV; and a = Atg eV/A and K = 32t 2 l cV/A 2 . We set 
the lattice constant to a = 1, and all the energies are in 
units of to. In real oxide materials, 13 ' to is in the range 
1.3-1.5. 

We also calculate the linear lattice dynamics from a 
renormalized electron Green's function within an adia- 
batic approximation. 12 ' From the eigenmodes of the lat- 
tice oscillation, we plot the corresponding dynamic neu- 
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tron scattering cross section: 14,15 ' 
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(a) HSD stripe state 
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(c) LSZ stripe state 
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(d) holes in HS and LS states 
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Fig. 1. (Color online) Different stripe states: (a) HSD, (b) LSD, 
and (c) LSZ. The circles (radius) and the arrows (length and di- 
rection) represent the hole and spin densities, respectively. Lines 
at the positions of O-ions represent static displacements. The 
shadows show the stripes. These are plotted for pressure param- 
eters in the range 1.33-1.35. The different hole states in the AF 
(in the HSD state) and NM (in the LSZ state) domains are de- 
picted in (d). The numbers represent the up- or down-spin hole 
densities calculated with A = 1.30 for HSD and A = 1.33 for 
LSZ. 



At ambient pressure (A = 1.0), the groundstate is 
calculated to be a diagonal stripe state, as experimen- 
tally observed. 8 ~ 10 ) The spin and charge configurations 
of the diagonal stripe state are shown in Fig. 1(a), where 
we observe the stripes with anti-ferromagnetic (AF) do- 
mains between them. Here, we term the state with such 
AF domains a high-spin (HS) stripe state in order to 
distinguish it from another stripe state which consists of 
almost non-magnetic (NM) domains [see Fig. 1 (b)]. We 
term the latter a low-spin (LS) stripe state. Details of 
these spin states will be discussed below. The stripes we 
find have diagonal, vertical, and zigzag characters. 

First, we show that the HS diagonal (HSD) and LS 
zigzag (LSZ) stripe states [Figs. 1(a) and (c)] are HF 
groundstates in our model. In Fig. 2, the energies of 
these states are plotted for different A, together with 
other stripe states: the LS diagonal (LSD), the HS verti- 
cal (HSV), and the LS vertical (LSV) stripe states. The 



periodicity of the stripes is 3 sites for all states; i.e., there 
is one hole per stripe. There is a transition from the HSD 
to the LSZ stripe state around A = 1.32. The LSD state 
has a close but slightly larger energy than the LSZ state. 

Next, we consider details of the transition. We first 
focus on the change of the stripe orientation and dis- 
cuss changes of the local spin state separately below. 
In our previous study on a C112O3 ladder system, 16 ' we 
found the tendency that the zigzag configuration be- 
comes preferable to the diagonal one under high pres- 
sure, although only the HS states appear in these cuprate 
systems. Therefore, we suppose that the transition from 
LSZ to LSD would not occur even when the pressure is 
increased further. In fact, the energy difference of these 
two states becomes larger as pressure increases and this is 
consistent with the behavior in the ladder system. How- 
ever, if a transition from HSD to LSD exists for some 
reason (e.g., different model parameters, different dop- 
ing, incomensurability, etc.), we would find the transition 
from LSZ to LSD at a higher pressure. 

Both the HSV and LSV stripe states cannot be the 
groundstate within the region we consider here. However, 
the LSV stripe state might become more stable in the 
higher pressure region by analogy with the CU2O3 ladder 
case. 16 ' Therefore, we discuss the properties of the HSD, 
LSZ and LSD stripe states, including the possibility of 
the HSD-LSD transition. 

There are distinct differences in the spin and hole 
densities between the HS and LS states. Clear differences 
are found in the AF or NM domain between the stripes 
[See Fig. 1(d)]. In the AF domain of the HSD state, one 
site, which has four local hole states (two d orbitals with 
two spins), is nearly half-filled: ~ 1.90 holes (« 48% fill- 
ing) at the site. The hole density on such a site is equally 
distributed into the two orbitals (d x 2_ y 2 and d^ 2 -^)- 
Here the hole states in the different orbitals on the same 
site are dominated by one spin: The dominant spin state 
has more than 97% of the holes in the orbital. The densi- 
ties and ratio between the local hole states do not change 
much for 1.0 < A < 1.32 

On the other hand, the NM sites in the LSZ state are 
occupied by ~ 1.80 holes each (ss 45% filling). Also, there 
is no longer equal distribution into the two orbitals: ~ 
85% in d x 2_ y 2 and ~15% in d^-r 2 - The up- and down- 
spin states are almost equally occupied: One of them has 
55% of the holes in the orbital. The LSD state has a 
similar hole distribution to the LSZ. 

As for the physical picture of the transition from 
HSD to LSD, it would be considered that the pressure 
enhances the hopping of holes and the system tends to be 
rather metallic. This effect should prefer the NM domain 
in order to gain the kinetic energy. The metal-insulator 
transition, however, cannot be treated in our model prop- 
erly due to the small system size. 

Thus far, we have shown that our calculation pre- 
dicts a transition between the HS (HSD) and LS (LSZ) 
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the path (0,0)-(27r, 0), the HSD case has additional modes 
around 61io meV; and the LSD case has high intensity 
in the range 55-65 to meV. 



1.25 1.3 1.35 

Pressure parameter \ 

Fig. 2. (Color online) Energies of different stripe states for dif- 
ferent A. The periodicity is 3 sites for all of the states: one hole 
per stripe. The lowest-energy state changes from HSD to LSZ at 
A ~ 1.32. 



states. Now we discuss how we can observe the transi- 
tion, taking into account the possibility of the HSD-LSD 
transition. X-ray diffraction should be able to detect a 
change in stripe orientation from diagonal to zigzag as a 
change of the charge order from (^p, ^L) to (f ,7r). 

However, it would be difficult to observe the change 
of the local spin density itself. Thus, if the transition 
does not change the stripe orientation (HSD-LSD), it is 
difficult to detect directly. Therefore, we consider the ob- 
servation of the transition from a different perspective, 
namely, lattice dynamics. Theoretical 14 ' 17 ^ and experi- 
mental 18,19 ) studies show that stripe states are accom- 
panied by characteristic phonons due to the inhomoge- 
neous structure of charge and spin densities coupled to 
ion positions. 

In Fig. 3, we plot the phonon spectra 5'(k, w)/|k| 2 
[Eq. (5)] of the HSD, LSZ, and LSD stripe states, cal- 
culated near the transition point. In the LS cases, the 
whole spectra are softened compared to those of the HSD 
state. The softening is ~ 7t meV (« 56to cm" 1 ) at (0, 0) 
and ~ 15t meV (« 121t cm" 1 ) at (tt, 0) in the LSD 
case. In the HSD case, phonon eigenmodes exist around 
61to, 65to, and 75-78 to meV near the transition point. 
At A = 1 (data not shown in figure), the correspond- 
ing branches to those at ~ 61to and 65 to meV are at 
~ 55to and 60to meV, respectively, and they are shifted 
upwards in frequency as A increases. The 75-78 to meV 
branch does not change much. 

In the LSD case, the phonon eigenmodes lie in the 
range 45-78 to meV. However, most of the modes have 
much weaker intensity than the strongest mode, which 
we will discuss later: In Fig. 3, we only see the spectra 
of high-intensity modes. 

Along the path (0,0)-(27r, 2tt), most of the intensity 
is around 65to and 75-78 to meV for the HSD case, and 
around 56to and 70-73 to meV for the LSD case. Along 
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Fig. 3. (Color online) Phonon spectra for HSD, LSZ, and LSD 
stripe states. Two different paths are plotted: (0,0)-(2,2) [(a), 
(b), (c)] and (0,0)-(2,0) [(d), (e), (f)], where n = 1. The results 
are symmetrized with respect to x and y directions. The pressure 
parameter is set as A = 1.30 for HSD [(a), (d)], 1.33 for LSZ [(b), 
(e)], and 1.345 for LSD stripes [(c), (f)]. The lines are guides to 
the eye. 



The highest intensity modes for the HSD and LSD 
cases are at 65to and 58to meV, respectively. Both modes 
are found at (0, 0) and show the same oscillation pat- 
tern [Figs. 4(a) and (b)]. These modes involve the O ions 
around the stripes and should be connected to the charge 
sliding of the diagonal stripes. (A similar mode has been 
calculated for cuprates in Ref. 14 - 1 — these are distinct 
from the breathing-type 18 -' modes.) 

The highest intensity mode in the LSZ is found at 
56tn meV, and this is also related to the stripe sliding 
[Fig. 4(c)]. In the LSZ case, however, the intensity of 
this mode is not so high as in the other cases, and there 
are several modes with similar intensity. For example, 
one of the modes involves all of the O y ions around the 
stripes, oscillating at 58to meV. There are some modes 
involving off-stripe ions, as well. We discuss the lower 
energy modes, other than the sliding one, below. 

We have seen that the high-intensity modes are the 
same sliding-related oscillation in both the HSD and LSD 
cases, and also found a similar mode in the LSZ. Here, we 
discuss the differences of the low-energy modes in the HS 
and LS cases, in the energy ranges 60-70 to meV (HSD) 
and 55-65 t meV (LSD, LSZ). 

In the HS case, all of these modes are spatially local- 
ized and involve only the O ions around the stripes. The 
vibration pattern of these modes is well-reproduced as a 
combination of three local vibrations (Fig. 5). The local 
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Fig. 4. (Color online) Lattice excitation pattern of high-intensity 
modes for (a) HSD case at ~65trj meV, (b) LSD case at ~58to 
meV, and (c) LSZ case at ~56fo meV. The filled circles represent 
the stripe sites. 



vibration patterns in Figs. 5(a) and (b) have the center- 
of-mass shift perpendicular and parallel to the stripes, 
respectively. The one in Fig. 5(c) is a Jahn- Teller (JT)- 
like vibration. 



(a) perpendicular (b) parallel 



(c) JT-like 
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Fig. 5. (Color online) Lattice vibration pattern around stripes. 
The thick line shows the stripe direction. The vibration pattern 
of local modes in the HS case are represented by a combination of 
these three patterns: (a) perpendicular to the stripe, (b) parallel 
to the stripe, (c) JT type. 



In the LS cases, the modes in the 55-65 to meV range 
include non-local modes as well, which involve all or most 
of the O ions. The phonons for the LS cases do not have 
as much of the original stripe character as those for the 
HS case, although there are still such modes, including 
the sliding one. It would be more accurate to say that the 
LS cases have some non-stripe modes in the lower-energy 
range, while the HS one has such modes but they are all 
in the higher-energy range, separated from the lower- 
energy stripe modes. However, since the correspondence 
of the modes is weak between the HSD and LSD states 
(except the sliding one), it is difficult to definitively as- 
sign each mode. 

From the analysis of the calculated phonon spec- 
tra, we conclude that the phonon excitations in the LS 
stripe state have less stripe signature than in the HS 
stripe state. From this, it can be deduced that the pres- 
sure reduces the stripe feature toward the inhomogeneous 
charge and spin structure in the high pressure limit. This 
is consistent with the physical picture of the HSD-LSZ 
transition discussed above. 

We now make an estimate of the transition point 



within a meanfield approximation. The hopping energy 
at the transition point corresponds to a lattice parameter 
change of about 8% (^ = A 33 « 1.08). From the Birch- 
Murnaghan equation of state 20 ' and the bulk modulus of 
NiO, 21 ' the corresponding pressure is estimated as about 
70-90 GPa. Although the estimated transition point is 
quite high, the actual value should be lower, since the 
estimate is within a meanfield approximation. 

In summary, we have investigated pressure effects on 
striped nickelate compounds by considering a four-band 
Hubbard model and allowing for a static displacement of 
oxygen ions. Using a selfconsistent Hartree-Fock calcula- 
tion, we found a pressure-induced transition from an HS 
stripe state to an LS one. The LS stripe state is char- 
acteristic to nickelate system, where two d orbitals are 
involved to make hole stripes, unlike those in cuprates, 
which involve only one d orbital. 

The prediction of this transition could be verified 
by observing a change in the x-ray diffraction pattern 
or in the phonon spectra. The former is useful if the 
state changes into zigzag stripes. If the stripe orienta- 
tion does not change at the transition, the latter can 
provide signatures of the transition. At the transition 
point, a softening of the lattice excitations around 60-65 
t meV (w 484-524 t cm" 1 ) by 7t meV (« 56t cm" 1 ) 
should be observable by means of infrared measurements 
under high pressure. Measurement by neutron scattering 
should also find a sharp change in the phonon spectra. 
However, the transition point may be too high for a neu- 
tron experiment. 

The phonon excitations in the LS stripe state have 
less stripe signature than in the HS stripe state. We spec- 
ulate that the LS stripe state may be metallic and that 
the transition has a metal-insulator transition character. 
If so, the optical conductivity should show a change near 
the transition point. 
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